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ABSTRACT 

Our 5 ks Chandra ACIS-S observation of the globular cluster NGC 6652 detected 7 X-ray sources, 
3 of which are previously unidentified. This cluster hosts a well-known bright low-mass X-ray binary, 
source A (or XB 1832-330). Source B shows unusual rapid flaring variability, with an average Lx(0.5- 
10 keV) ~ 2 X 10"^^ ergs/s, but with minutes- long flares up to Lx = 9 x 10^"* ergs/s. Its spectrum can 
be fit by an absorbed power-law of photon index F ^ 1.24, and hardens as the countrate decreases. 
This suggests that part or all of the variation might be due to obscuration by the rim of a highly 
inclined accretion disk. Sources C and D, with Lx ~ 10'^^ ergs/s, have soft and unusual spectra. 
Source C requires a very soft component, with a spectrum peaking at 0.5 keV, which might be the 
hot polar cap of a magnetically accreting polar cataclysmic variable. Source D shows a soft spectrum 
(fit by a power-law of photon index ~ 2.3) with marginal evidence for an emission line around 1 keV; 
its nature is unclear. The faint new sources E, F, and G have luminosities of 1-2x10'^^ ergs/s, if 
associated with the cluster (which is likely). E and F have relatively hard spectra (consistent with 
power-laws with photon index ~'1.5). G lacks soft photons, suggesting absorption with Nh > 10^^ 

_2 

cm . 
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1. INTRODUCTION 

The high-density environments of globular clusters 
were suggested to be caus es of X-ray binary pr oduction 
early in X-ray astronomy (|Katzlll975t IClark|[l975iy Sev- 
eral types of X-ray source s have now been identif ied in 
Galactic globular clusters (jVerbunt fc Lewinll2004[) . Fif- 
teen luminous low-mass X-ray binaries (LMXBs) reach 
Lx > lO'^^ ergs/ s, some only during short outbursts 
(iSidoli et al.|[200il: IVerbunt fc 'Lewi]3[200l iHehlke et al.l 
120101: iPoolev et all 1201011 The low-luminosity X-ray 
sources Lx = 10^^ - lO^'' ergs/s include several types 
of systems such as quiescent LMXBs (containing neu- 
tron stars between accretion episodes), millisecond ra- 
dio pulsars, cataclys mic variables (CVs), an d magneti- 
cally active binaries (jVerbunt fc Lewinll2004[ ). Quiescent 
LMXBs containing neutron stars (qLMXBs) often show 
soft t hermal spectra dominated by blackbody-l ike emis- 
sion (iRutledge et al.ll2002al: [Heinke et al.ll2003D . though 
a spectrally harder component is often present and 
sometimes dominant (ICampana et al.ll200^lJonker et al.l 
I2004bt IWiinands et al.ll2005D . The origin of the soft com- 
ponent is generally thought to be thermal X-ray emission 
from the heated NS surface, modified by the star's hydro - 
gen atmosphere (jZavlin et al.lll996t iBrown et al.l 119981 ). 
The harder power-law component is not yet well under- 
stood, but may be prod uced by continuing ac cretion or 
possibly pulsar activity (jCampana et al.|[l998[ ). 

The globular cluster NGC 6652 c ontains one luminous 
{Lx > 10^6) LMXB, XB 1832-330 (jPredehl et al.l[T99l 
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and sever al lower-luminosity X-ray sources. This core- 
coll apsed (iNovola fc Gebhardtll2006f ) cluster is ~11.7 Gyr 
old (jChabover et al.ll2000[) . 9."0±0.4 kpc from the sun and 
suffers extinction with equival ent hydrogen column den- 
sity Nh = 5.0 X 10^0 cm-2 (Harris 19961). The three 
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fainter sources detected by iHcinkc et al. 
ks Chandra HRC-I observation are not well understood. 
The brightest of the three low-Lx sources (source B) was 
seen at luminosities of a few 10'^^ ergs/s while its opti- 
cal counterpart lies on the main sequence (1 magnitude 
below the turnoff) in aV — I optical color-magnitude d i- 
agram, suggesting a qLMXB nature (jHeinke et al.l l200ll . 
However, the optical counter part has also b e en ob served 
to show blue U-V colors (jPeutsch et al.l I1998D . and 
strong variability wi th a (possible) 43.6-minute period 
(jPeutsch et al.l l2000f ). which, if it is the true period, 
would exclude a main-sequence companion. Source C 
has a very blue optical counterpart, which indicates a 
bright disk and a relatively high rate of mass transfer. 
Combining this with its relatively low Lx, the efficiency 
of energy extraction is inferred to b e relatively low, sug - 
gesting a CV rather than an LMXB (jHeinke et al.ll2001| i. 

In this paper we present a 2008 Chandra ACIS-S ob- 
servation of NGC 6652's X-ray sources. We also consult 
2000 Chandra HRC-I data and 1994 ROSAT HRI data 
for long-term variability information. 

2. DATA REDUCTION 

We observed NGC 6652 on June 9, 2008 for 5.6 ks 
with Chandra 's ACIS-S detector in a 1/4 subarray 
mode. The data were reduced using the CIAO version 
4.3 softwarcQ. We created a new bad pixel file with the 
acis_runJiotpix script. The level 1 event files were re- 
processed by calibrating for charge-transfer inefhciency 
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on the detector and time-dependent gain adjustments. 
The data was then filtered for grade and status bits ac- 
cording to the standard CIAO Science Threadfl A back- 
ground hght curve shows no evidence of flaring. 

The CIAO WAVDETECT program was run using the 
energy range 0.3-7.0 keV. We find 7 sources, all of which 
are locat ed within (or at) the cluster half-mass radius 
rh, 0.'65 (|Harrij[T996[) . Three faint sources (E, F and 
G) were not visible in the HRC observation of NGC 
6652 by iHeinke et al.l ()2001[ ). The positions of sources 
B through G, shown in Table 1, were computed using 
WAVDETECT centroiding. Due to the high flux of pho- 
tons from source A, pileup effects produce a characteris- 
tic "donut hole" at the source location and a prominent 
readout streak, which prevents WAVDETECT from ac- 
curately computing the position of source A. We there- 
fore estimated the location of source A by matching a 
symmetric circle to the "donut hole" , by eye. Running 
WAVDETECT in the energy bands 0.3-2.0 keV, 2.0-5.0 
keV, and 5.0-7.0 keV did not locate additional sources. 
The positions of sources A, B, C and D are consistent 
with those from the 2001 Chandra observation. We show 
a smoothed ACIS image of N GC 6652 in Fig 1. U sing 
the log N-log S relationship of iGiacconi et al.l (|2001| ) , we 
calculated the expected number of active galactic nuclei 
observed within one half-mass radius of the cluster to be 
0.14, indicating that all detected X-ray sources are likely 
to be members of the cluster. 

Due to pileup effects, the data from the source re- 
gion of A is not reliable. Pileup occurs when two or 
more photons are recorded by the CCD as a single 
higher energy event, which results in distortion of the 
obser ved source s pectrum and an underestimated count 
rate (|Davislf2001[ ). The high count rate also produces a 
readout streak, as photons land on the detector during 
the short frame readout time. By extracting a spectrum 
from the readout streak, we are able to obtain spec- 
tral data from source A without having to model the 
pileup. Using dmextract we extract a source spectrum 
from the readout streak, and a background spectrum 
from surrounding source-free regions of the observation. 
Wc then create new RMF and ARF files using the tools 
mkacisrmf and mkarf respectively, and correct the streak 
exposure time to 110.4 s to account for the fact that the 
spectrum was extracted from the ACIS readout streak. 

For sources B, C and D, we extracted 0-10 keV spectra 
using specextract from circular regions, of radius 2" for 
B and 1.4" for C and D, located at the source positions 
produced by WAVDETECT. Wc produced spectra, re- 
sponse files and backgrounds following the appropriate 
CIAO threacQ, correcting the response files for the frac- 
tion of the flux contained within the extraction aperture, 
group the spectra, and fit them within XSPEC0. Count 
rates were calculated for sources E, F and G directly 
from the level 2 event file using the dmstat tool, since 
these sources are too faint to allow spectral analysis. The 
count rates were then converted into unabsorbed fluxes 
with PIMM^, by assuming that the faint source spec- 
tra can be modeled by a power law of photon index F 

® http:/ /cxc. harvard.edu/ciao/threads/all. html 

http:/ /cxc. harvard.edu/ciao/throads/psextract/ 
* http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xspec/ 
^ http:/ /asc. harvard.edu/toolkit/pimms.jsp 
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Fig. 1.— Chandra ACIS-S image of NGC 6652. The data is 
smoothed with a gaussian kernel radius of 1 " . The 7 detected 
sources are labeled from A to G. Both the core radius (0.07' ) and 
the half mass radius (0.65' ) are indicated. 



= 1.4, typical for globular cluster sources at this Lx 
(|Heinke et al.l [20051 ) . and consistent with their hardness 
ratios (below). 

To study long term variability of the sources, we re- 
analyze the archival HRC-I data from a 2000 observa- 
tion of NGC 6652, filtering the data for grade and sta- 
tus. We extract hghtcurves for sources A, B, C, and D 
with the CIAO tool dmextract and estimate source lumi- 
nosities with PIMMS. Wc extract counts and lightcurvcs 
from source regions of radius 2" centered at the WAVDE- 
TECT coordinates obtained from our ACIS observation. 

3. X-RAY ANALYSIS 

3.1. Source A 

The 462 events extracted from source A's readout 
streak were grouped into 9 bins, with 50 counts per 
bin to maximize spectral resolution while minimizing the 
error in the normalized counts in each bin (Fig. 2). 
Wc ignore events above 8 keV, as the ACIS response 
is poorly understood at these energies. The spectrum 
is well fit by an absorbed power-law model with a pho- 
ton index of F = 1.7^!^-^ (Table 2). The Nh value of 
2.7l°:l X 1021 cm exceeds the accepted cluster value 
of 5x10 cm' -2 (lHarrislll996n . We der ivc an unabsorbed 
luminosi ty for source A of ( 0.5-6.0) = 4.4t°-^ x lO^^ 
crgs/s. iMukai fc Simli ()2OOO0 observed XB 1832-330 
with ASCA, finding that a partial covering model was 
required, instead of a single absorption column, to fit the 
spectrum. A similar model was r equired by Bcppo SAX 
(jParmar et al.l[200l and XMM (jSidoh et al.l[2008h ob- 
servations. Our spectrum is of insufficient quality to dis- 
tinguish between a partial covering model vs. a single 
absorption column. 

We detect 9205 counts for source A in the 2000 HRC 
observation of NGC 6652. We convert the average source 
count rate into an unabsorbed bolometric X-ray luminos- 
ity with PIMMS using our best ACIS spectral fit, finding 
Ljf (0.5-6.0) = (1.61 ± 0.02) xlO^^ ergs/s. A's luminosity 
appears to ha ve decreased by alm ost a factor of four from 
2000 to 2008. iSidoh et al.l H2008D report two XMM obser- 
vations in Sept. and Oct. 2006. Converting their X-ray 
luminosities to the 0.5-6 keV range, they find Lx{0.5- 
6 keV )=1.2 X lO^s and 1.1 x lO^^ergs/s. iParmar et all 
(|2001f) report a BeppoSAX measurement in March 2001, 
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Fig. 2. — Top: Chandra 2008 X-ray spectrum of source A con- 
taining 9 bins of 50 counts per bin. The spectrum is best-fit by an 
absorbed power- law model of photon index of F = l.T^g'j- Bot- 
tom: Residuals to the best fit. 

which we convert t o Lx(0.5-6 keV)=1.5 x 10"^^ ergs/s. 
iTarana at all ()2007D report 2003-2005 INTEGRAL mea- 
surements (restricted to >20 keV, and thus not directly 
comparable to our fluxes) that are consistent with the 
BeppoSAX flux. RXTE PGA Galactic Bulge Scan mon- 
itoring of XI832-33Cl3 indicates substantial variability 
on ~6-montli time scales, but also a general declining 
trend from 1999 to 2010, by about 30%. 

3.2. Source B 

3.2.1. Timing Analysis 

We produced a lightcurve by binning the 0.3-7 keV 
data into 103 bins each with 50s per bin (Fig. 3), which 
shows clear variability, by factors >10 on timcscalcs <100 
s. (A K-S test indicates variability at >99.9% confi- 
dence.) Such variability is unusual in low- mass X-ray 
binaries. To test whether the variation may be caused 
by changes in obscuring column (perhaps, if the system 
is edge-on, caused by material at the rim of an accretion 
disk) , we group the lightcurve by countrate, with bound- 
aries when the countrate crosses 0.15 and 0.3 counts/s 
(as the count statistics are low for 50 s bins). Within 
each larger bin, we compute the hardness ratio of 0.3-1 
keV photons (the most likely to be absorbed) over 0.3-7 
keV photons (Fi g. 3, bot t om), with binomial la error- 
bars derived from lGehreld ()1986f ). If the dips were due to 
obscuration, we might expect the lowest-countrate bins 
to have the lowest hardness ratios. This is suggested by 
the first 1000 s, but is less clear from the dataset as a 
whole. We quantify this by comparing the ratio of 0.3-1 
keV counts to total counts in the low-countrate portions 
(<0.15 cts/s), vs. medium and high-countrate portions 
(below and above 0.3 cts/s), finding the fraction to be 
0.15±0.03 for the low-countrate portions vs. 0.25±0.04 
and 0.23±0.03 for the higher-countrate portions. This is 
a significant effect, indicating that obscuration may play 
a role in the dipping. 

The peaks in the lightcurve reach countrates that sug- 
gest substantial pileup. We estimate the amount of 
pileup and intrinsic luminosities using the PIMMS tool, 
for the best-fit absorbed power-law spectrum (see be- 

http: / /Iheawww.gsfc. nasa.gov/users/craigm/galscan/html/R_1832- 
330.html 



low; the choice of spectrum does not appreciably affect 
results). At the peak observed countrate of 0.46 cts/s, 
the estimated fraction of recorded events that are actu- 
ally multiple is 27%, and we infer the intrinsic Lx(0.3- 
7)~ 9 X lO'^^ crgs/s, a factor of 2 higher than would be 
extrapolated without pileup. However, the pileup model 
in PIMMS is systematically uncertain, so the true peak 
Lx is uncertain, presumably by less than a factor of 2. 
The inferred minimum is < 2 X 10^3 ergs/s, for 0.02 
cts/s (cf. Fig. 3). 
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Fig. 3. — Top: ACIS countrate lightcurve of source B, 0.3-7 keV, 
in 50 s bins. Lx/ctrate conversion estimated as 1.2 X 10^* ergs/s 
for 0.1 cts/s, reaching 9 X 10'^'' ergs/s for 0.46 cts/s, due to pileup. 
Bottom: Ratio of counts in 0.3-1 keV vs. those in 0.3-7 keV. 
Boundaries of ratio bins set by when countrate crosses 0.15 and 
0.3 keV. Bins over 0.3 cts/s: (red) filled squares; 0.3-0.15 cts/s: 
(magenta) open triangles; and under 0.15 cts/s: (blue) crosses. 

Power spectra (produced using 

XRONO^B) show no 
evidence of periodicity during the 5.6 ks AGIS obser- 
vation. We re-examined the 2000 HRG observation to 
check for long-term variability. Using PIMMS and our 
best AGIS spectral fit (below), we estimate the unab- 
sorbcd luminosity to be Lx(0. 5-6.0) ~ (1.1 ± 0.1) x 
10^4 ergs/s, vs. our average AGIS unabsorbed luminos- 
ity of Lx(0.5 - 6.0) = X 10^4 ergs/s (section 
3.2.2). iHeinke et al.l (j2001l ) found significant variability 
during the HRG observation. The short, low-sensitivity 
lightcurve suggests variation by at least a factor of two. 
As the HRG sensitivity and exposure length are each 3 
times smaller than those for the AGIS observation, its 
usefulness for studying variability is limited. 

An archival ROSAT HRI observation (March 27, 1994, 
for 817 s exposure time) shows marginal evidence for 
source B. The image (Fig. 4) shows one bright source, 
which we attribute to A. setting our astromctry. We mea- 

http:/ /heasarc. nasa.gov/docs/xanadu/xronos/xronos. html 
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sure 6 photons within 2" of source B's position, while 
only lA^I'^ photons are expected at this position, us- 
ing the average background rate between 6.5" and 10" 
from source A, giving 4.6^2.5 (^^ ^ sigma; 4.6^4^5 
2 sigma, iGehreljfTOSl) photons from B. This ~2cr de- 
tection suggests (correcting for the HRI point-spread 
function enclosed energy of 38% within 2") (0.5-6 
kcV)= 1.5lg g X 10'^'* crgs/s. Thus, we have no evidence 
for B's variability on long time scales. 




Fig. 4.— ROSAT HRI image, March 1994, of NGC 6652. The 
(corrected, using A's position) positions of sources A, B, C and 
D are marked with small (radius=2") circles, while the annulus 
around A used to estimate the background is indicated with larger 
circles (radius 6.5" and 10"). 

3.2.2. Spectral Analysis 

We created both an overall time-averaged spectrum of 
B (680 counts, binned by 30 counts/bin), and spectra 
for the three countrate ranges lis ted above. W e include 
the Chandra CCD pileup model (|Davi!j|200lD . with the 
grade morphing parameter a fixed to 0.5, and photoelec- 
tric absorption (XSPEC phabs). B's time-averaged spec- 
trum can be fit well with an absorbed power-law model 
with photon index T = 1.24^°;^° (Table 2). We also try 
adding a NS hydro gen atmosphere model (NSATMOS, 
iHeinke et al.l l2006a[ ) to the absorbed power-law model. 
We fix the mass and radius to canonical values (1.4 Mq, 
10 km), so kT is the only free parameter. The NSATMOS 
component does not improve the quality of the spectral 
fit, and the upper limit on its kT implies an upper limit 
on Lx,Nsi0.5-l0.0 kcV) < 7x 10^^ ergs/s. This limit 
may not be applicable if the system is edge-on, as the 
NS could suffer higher Nh than other X-ray emitting 
regions. 

We next fitted the three countrate-selected spectra, 
with 240, 186, and 253 counts [low to high], binned by 
15 counts each, to investigate what parameters may be 
changing. An absorbed, piled-up power-law spectrum 
is best fit with (at least) two parameters varying; the 
power-law normalization, and either the Nh or power- 
law spectral index. We give parameters for both such fits 
in Table 3. The alternative of allowing only the power- 
law normalization to vary gives significant residuals in 
the spectra (Fig. 5), and a higher chi-squared: 42.93 for 




Energy (keV) 

Fig. 5. — Top: Chandra X-ray spectra of source B at high (black, 
crosses), medium (red, diamonds) and low (green, filled squares) 
countrates. All are fit with a power-law model with only normal- 
ization varying (see text). Bottom: Residuals to this fit, demon- 
strating the spectral changes. 



38 degrees of freedom (dof), vs. 30.3 or 30.6 (for 36 dof) 
from the other two models (respectively) . An F-test gives 
probabilities of 0.2% of attaining such an improvement in 
the chi-squared by chance, indicating that another quan- 
tity besides the normalization is also varying. 

3.3. Source C: 

An unbinned 0.3-7 keV lightcurve gives a 0.01 K-S 
probability of source constancy, indicating that C is vari- 
able. This is confirmed by its (0.3-7 keV, 500 s binning) 
lightcurve, which shows two clear peaks (Fig. 6). A 
power spectrum shows no evidence of periodicity. Us- 
ing the best- fit double mekal spectral model (below), 
we find a peak Lx(0. 5-10.0) = 3.41;^;} x lO^^ ergs/s, 
and a minimum Lx < 3.2x 10"^^ ergs/s, for an average 
Lx = 1.1 X lO'^"^. Using the power- law spectral model (a 
poor fit) to the ACIS data to infer the HRC spectrum, 
wc find Lx = 2 X 10^^ crgs/s, half the ACIS estimate 
with this spectrum, but the spectral uncertainties render 
this conclusion uncertain. 

We group C's spectrum by 10 counts/bin for statis- 
tics (Fig. 7), and also fit the unbinned spectrum with 
C-statistics (finding similar results). No single model we 
tried can fit this very soft spectrum, with reduced 
values well above 2 (Table 2). The spectrum can be fit 
by a model contai ning two MEKAL (thermal plasma, 
iLiedahl et al.lll995() components with the cluster absorp- 
tion (to simplify fitting we fixed Nh), with temperatures 
of < 0.096 and > 2.3 keV. Ahhough a double MEKAL 
model is o ften used to describe the spectra of CVs (e.g. 
iBaskill eTal. 2005) . such a strong low-temperature com- 
ponent is rarely seen. We are aware o f such components 
only in nova remnants ()Balmanll2005D . Alternatively, an 
absorbed blackbody and MEKAL repres ents a simpli- 
fied f orm of a model for polar CVs (e.g. iRamsav et al.l 
l2004f l. C's spectrum is reasonably fit by a blackbody 
plus MEKAL model, with a blackbody temperature of 
Q7l:lt cV and inferred radius of 46t^^ km, and a MEKAL 
temperature of >3.2 keV. 

3.4. Source D: 

A K-S test on the unbinned lightcurve gives a prob- 
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Fig. 6. — Chandra lightcurve of source C, 0.3-7 keV, at 500 s 
binniiiEf. 
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Fig. 7. — Top: Chandra 2008 X-ray spectrum of source C, fit by a 
low-temperature blackbody plus MEKAL model. Dotted lines in- 
dicate the two components (the blackbody is the lower-temperature 
component). Bottom: Residuals to the best fit. 

ability of constancy of 0.46. The estimated luminosity 
from the HRC data using the ACIS power-law fit below 
is Lx(0. 5-10.0) = l.lx 10^^ ergs/s. This is comparable 
to the luminosity from the ACIS data, Lx (0.5-10.0) = 



i+5 



10 ergs/s, so no variability is apparent. 



We extract a 69-count spectrum, binned with 10 
counts/bin. Source D is relatively soft, with all 6 bins 
below 2.0 keV (Fig. 8), motivating us to freeze Nh to 
the cluster value. An NSATMOS model is a poor fit, but 
an absorbed powerlaw (photon index of 2.3l;Q g) is an ad- 
equate fit to the binned spectrum; a single-temperature 
MEKAL is a somewhat worse fit (Table 2; null hypothesis 
probabilities of 25% and 11% respectively). The power- 
law fit, however, shows residuals suggesting strong line 
emission around 1 keV (Fig. 8). Adding a gaussian 
line of zero width to the power-law fit gives a line en- 
ergy of 1.03lQj!)4 keV, but does not improve the fit by 
a statistically significant amount (an F-test indicates a 
probability of 43% that such an improvement could hap- 
pen by chance). Fitting the unbinned data (using the 



0.2 0.5 1 2 

Energy (keV) 

Fig. 8. — Top: Chandra spectrum of source D, fit to an absorbed 
power-law model (see text). Bottom: Residuals to the best fit. 

C-statistic) with a power-law gives a slightly larger pho- 
ton index (2.7^0 4), and the fraction of simulated spectra 
with lower C-statistic values is 97%, suggesting a rel- 
atively poor fit. (The single-temp MEKAL fit to the 
unbinned spectrum gives 96% of simulations with lower 
C-statistic values, as well.) 

3.5. Faint Sources: 

We detect only 15, 7, and 15 events for source E, F, and 
G respectively, too few for detailed spectral and timing 
analyses. We extract background-subtracted count rates 
in energy bands of 0.5-6.0 keV, 0.5-1.5 keV, and 1.5-6 
keV, using the WAVDETECT regions and PIMMS to 
estimate unabsorbed source luminosities. While sources 
E and F have roughly equal numbers of soft and hard 
counts, all but one of source G's counts lie above 1.5 
keV, suggesting strong absorption. 

Using PIMMs, we estimate the unabsorbed luminosi- 
ties of the faint sources in both the 0.5-1.5 keV and 1.5- 
6.0 keV energy bands, using an absorbed power-law of 
photon index 1.4. We determine the total unabsorbed 
luminosity from 0.5 to 6.0 keV for each source by adding 
the soft and hard luminosities (Table 1). We extract the 
counts at the ACIS WAVDETECT coordinates for each 
source from the HRC dataset, finding 4, 1, and counts 
from E, F, and G respectively. These upper limits are 
consistent with their measured ACIS luminosities. 

3.6. X-ray Color-Magnitude Diagram: 

We create an X-ray color-magnitude diagram (Fig. 
9) containing the sources observed in the 2008 ACIS 
data, by plotting the unabsorbed 0.5-6.0 keV luminosi- 
ties versus the X-ray color, defin ed as 2.5 log [(0.5-1. 5 
kev counts)/(l. 5-6.0 keV counts)] (iGrindlav et al.ll2001l) . 
We plot an X-ray color-luminosity relation for the NSAT- 
MOS model, and X-ray color predictions for power-law 
and MEKAL spectral models. The spectra of qLMXBs 
are frequently m odeled with a neutr on star atmosphere 
(e.g., NSATMOS fHeinke et alll2006b[ ) component, with a 
harder power-la-w comp onent with a photon index of 1-2 
(jRutledge et al.l[2"002b[ ). makin g up anywh ere from < 10 
to > 90% of the 0.5-10 keV flux (| Jonker eTa l. 2004aj). We 
plot a representative NSATMOS plus power-law model in 
which the power-law component contributes 50% of the 
0.5-6.0 keV luminosity. Note that B's position is shifted 
to the left by pileup (considering its rapid variability, it 
is hard to correct for this). 
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Fig. 9. — X-ray color-magnitude diagram for all sources detected 
in the 2008 Chandra ACIS-S observation. The color is defined as 
a function of the ratio of low-energy counts to high-energy counts. 
X-ray luminosities between 0.5 keV and 6.0 keV are plotted versus 
the color along with their respective errors. Also shown are the 
theoretical cooling tracks for the power-law, MEKAL, NSATMOS, 
and NSATMOS-l-power-law models. The NSATMOS-|-power-law 
model is defined such that 50% of the model's 0.5-6 keV luminosity 
is produced by the power-law component. 

No sources lie near the NS atmosphere coohng track, 
but sources C and D have colors and lum inosities in the 
range of qLMXBs in oth er clusters (e.g. iPoolev fc HutI 
l2006HHeinke et al.ll2006"bl ). C's rapid variability and un- 
usual spectrum likely rule out a qLMXB nature for it, 
but D could be a qLMXB with a dominant power-law 
spectral component. 

Surveys of X-ray sources in globular clusters indicate 
that most nonmagnetic CVs display hard spectra, typi- 
cally consiste nt with MEKAL temperatures > 6 keV o r 
colors < 0.5 (|Grindlav et al.llMlt Foolev et all [20021 ) . 
Sources E, F and G are consistent with the X-ray lumi- 
nosities and sp ectra of CVs observed in many other glob - 
ular clusters (jPoolev fc HutI 120061: iHeinke et ai]|2005[ ). 
though quiescent LMXBs, active binaries, and millisec- 
ond pulsars cannot be ruled out. G's extremely hard 
spectrum, though based on few counts, indicates a high 
intrinsic absorption (> 10^^ cm~^), and thus an edge- 
on CV (e.g. W8 W15, W33, and AK09 in 47 Tuc, 
IHeinke et al1l2005[) or a background AGN. 

4. DISCUSSION 

4.1. Source B: Unusual LMXB 

Source B, one of the brightest \ow-Lx cluster sources, 
can be classified as a very faint X-ray transient ( VFXT) . 
VFXTs are X-ray transients that have peak luminosi- 
ties of 10^"'"^^ ergs/s and quiesce nt luminosities at least 
one order of mag nitude lower (jWiinands et al.l 120061 : 
iMuno et al.l l2005aD . generally containing an accreting 
neutron star or black hole. Some VFXTs have "nor- 
mal" outbursts as well as very faint outbursts, but it is 



not cle ar if all do, or if VFXT beha vior has a variety of 
causes (jPegenaar fc Wiinandsl l20lol ) . 

Fig. 3 illustrates that source B reaches peak X-ray lu- 
minosities up to La-(0. 5-10.0 keV) = 9 x lO'^^ ergs/s and 
minimum Lx < 2 x 10'^^ ergs/s on timescales of min- 
utes. B's high peak luminosity is strong evidence that 
the system must contain a neutron star or a black hole, 
but the variability is unusual for LMXBs. Perhaps the 
simplest explanation for this variability is a high inclina- 
tion angle, so that we observe B's accretion disk edge-on, 
and our view of the central X-ray source is interrupted 
by structures at the ac cretion disk rim (|White fc Holtl 
119821 : iXiang et al.l l2009| ) . Obscuration by an accretion 
disk should lead to changes in the Nh value, and our 
spectral analysis gives evidence in favor of this (section 
3.2.2). Our spectral fitting requires intrinsic Lx changes 
along with Nh changes, which seems to argue against 
this explanation. However, our spectral fitting is con- 
strained by a lack of data to only three countrate ranges, 
which contain substantial variability, and we think it 
likely that variability also occurs on timescales shorter 
than we can accurately probe. These two effects could 
prevent us from accurately measuring how Nh changes 
with Lx in B's lightcurve, and thus we cannot yet rule 
out that obscuration is fully responsible for the variabil- 
ity 

If we only see scattered light from a central source (in 
this picture the variability is due to obscuration of an ac- 
cretion disk corona), the true isotropic X-ray luminosity 
could be higher t han the observed lu minosity, by a factor 
of up to 100 fc.g. lMuno et al.ll2005bf ). Thus, B could be a 
'normal' LMXB, which would remove the difficulty of ex- 
plaining its unusually low accretion luminosity. We have 
been awarded a deeper Chandra observation of NGC 
6652 (50 ks, to be taken mid-2011), which should clar- 
ify the spectral variability of this source. We have also 
been awarded a Gemini time-series imaging observation 
of NGC 6652, designed to search for evidence of short 
periods i n sources A an d B suggested by previous HST 
imaging (jPeutsch et al.l l2000: Hcinkc et al. 2001). 

4.2. Source C: A Second Polar in a Globular 
Cluster? 

The very soft spectral component of C is unlike any- 
thing seen in quiescent LMXBs, but is consistent with 
the spectra of some CVs. The (low-T) blackbody -I- 
(high-T) MEKAL spectral fit suggests a polar CV na- 
ture, where the blackbody describes soft X-rays emit- 
ted from the white dwarf's polar cap, and the MEKAL 
descri bes hard X-rays from the accretion column shock 
front (jRamsav fc Cropperll2004f ) . The inferred polar cap 
radius of 251^° km (from the blackbody spectral fit) com- 
pares well with the typical radius of an accreting pole o f 
a WD of -75 km (for a 0.6 MgWD, llshida et al.l[T997l) . 
The spectrum and luminosity are similar to the likely po- 
lar C V XIO in the globular cluster 47 Tuc (jHeinke et al.l 
|2005[ ). which has Ljf (0.5-6 keV)= 2.6 x lO^^ ergs/s, and 
was modeled with a kT = 53 eV blackbody and two 
MEKAL plasmas of temperatures 0.39 and >14 keV. 
Thus we suggest source C may be the second polar CV 
identified in a globular cluster. Knowing the frequency of 
magnetic accretion channelling in cluster CVs (by identi- 
fying magnetic CVs) is im portant for understanding their 
unusually rare outbursts ()Grindlavlll999t iDobrotka et al.l 
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[200I Hwova et al.|[2006l ). 

Alternatively, a very low-T thermal plasma has been 
seen in historical novae, from the nova shell (e.g. iBalmanI 
I2005| ). However, the nova eruption would probably have 
been seen if recent, or if old, the shell should have ex- 
panded to a resolvable size (e.g., the Nova Per 1901 shell 
would be 2.5" in radius if located in NGC 6652). Our up- 
coming Chandra observation should provide sufficiently 
high-quality spectra to distinguish between the spectral 
models discussed, and search for spectral variations over 
time. 

4.3. Source D: unusual CV or quiescent LMXB? 

Source D's luminosity and spectrum are consistent 
with both CVs or quiescent LMXBs. The feature of par- 
ticular interest is the hint of an emission line around 1 
keV. If this feature is real (which will be tested by our up- 
coming Chandra observation) , it would make D a rather 



odd X-ray source. There is another X-ray bright CV with 
a similarly strong low-energy line, the luminous CV X9 
in 47 Tuc (iHeinke et al.ll2005D . 
4.4. Faint Sources 

The pronounced spectral hardness of G indicates it is 
an edge-on CV, or a background AGN. E and F cannot be 
clearly classified, but as the majority of globular cluster 
X-ray sources with t heir X-ray colors and Lx are CVs 
(jPoolev fc Hutll2006[ ). we may suggest that CVs are the 
most likely possibility. Deeper Chandra observations 
should detect additional sources, enabling a comparison 
of the amazingly rich X-ray population of this globular 
cluster with other globular clusters. 

We acknowledge support from NASA Chandra grants, 
an NSERC Discovery Grant, and an Alberta Ingenuity 
New Faculty Award. We thank M. Muno for his assis- 
tance in obtaining these observations. 
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TABLE 1 
Sources in NGC 6652 Field 



Name 


Source 

CXOGLB J 


R.A. 


Err 


Decl. 


Err 


Counts 
(0-10 keV) 


Lx 

(0.5-1.5 keV) 


(10^2 ergs s 
(1.5-6 keV) 


(0.5-6 keV) 


A 


183543.7-325927 


18:35:43.69 


0.02 


-32:59:26.5 


0.5 


462" 




2813i|ig 


4386ltg«7 


B 


183544.6-325938 


18:35:44.567 


0.001 


-32:59:38.36 


0.02 


680 


271^ 


77 ± 8 


105+1° 


C 


183545.8-325923 


18:35:45.755 


0.003 


-32:59:23.21 


0.04 


90 


3.9lg:^ 


1.9+?« 


5 8+2-8 


D 


183545.6-325926 


18:35:45.648 


0.004 


-32:59:26.10 


0.05 


73 


q O + 0.8 


^■^-0.9 


5 4+1-8 


E 


183546.2-325929 


18:35:46.206 


0.010 


-32:59:29.26 


0.09 


15 


n 47+0-28 
'J-*' -0.19 


1 9+0.8 


1 6+1° 
l-"-0.7 


F 


183547.8-325920 


18:35:47.777 


0.004 


-32:59:19.84 


0.15 


7 




u.dO_Q 32 


80+°-*° 


G 


183548.8-325921 


18:35:48.816 


0.009 


-32:59:21.13 


0.10 


15 




r, 7+0.9 
^- '-0.7 


^•'-0.8 



Note. — Positions (with relative errors, not including systematic uncertainties on the Chandra astromctry), total counts, and inferred 
luminosities for NGC 6652 sources. 
From readout streak. 



TABLE 2 
Spectral Fits 



Src 


Model 


N„ 
(102° cm-2) 


r 


kT 

keV 


xS/dof 


Lx 


A 


POW 


27+1° 


1 7+0.3 
-'-•'-0.2 




0.87/6 


5.8+0 7 X 10^^ 


B 


POW 


c: + 6 
0-0 


^■•^-0.2 




0.70/18 


1.7+°-^ X 10^4 


B 


POW+NSATMOS 


5+8 
0-0 


1 9+0.2 
-■-•^-0.2 




0.736/17 


1.7l°:2 X 1034 




NSATMOS 






10+1°° 




< 7 X 10^2 


C 
C 


POW 
MEKAL 


(5) 
(5) 


5.2 


0.14 


3.14/5 
6.24/5 


4.1x103^ 
5.6x1032 


c 


MEKAL+MEKAL 
2nd MEKAL 


(5) 




0.08l°-°i5 
> 2 


0.45/3 


l.ltli X 1033 
7.3i4;3 X 1032 


c 


MEKAL+BBODY 


(5) 




> 3.3 


1.00/3 


1.2l:o;i X 1033 




BBODY 






n 067+°-°l3 




3.lj:i;3 X 1032 


D 


POW 


(5) 


9 .1+U.B 
^••-'-0.6 




1.33/4 


8+3 X 1032 


D 


MEKAL 


(5) 




9 + 5 
^-1 


1.88/4 


6.8l^;2 X 1032 


D 


NSATMOS 


(5) 




0.1 


2.3/5 


4.3 X 1032 


D 


POW+GAU 
GAU 


(5) 


9 q+0.7 
^••-'-0.7 

1 03+°i2 
^■^■^-0.04 




1.14/2 


8+1 X 1032 

5+1 X 1031 



Note. — Spectral fits to NGC 6652 sources (see text for details). Errors are 90% confidence for a 
single parameter, and are not computed if xt > 2. Luminosities are unabsorbed in erg s^i for 0.5-10 keV. 
Two-component models are continued on a second line (omitting the source name in the first column) , with 
the total luminosity on the first line and the second component's luminosity individually on the second 
line. All models include PHABS; for faint, soft sources, we have fixed the Nh to the cluster value. 
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TABLE 3 

Countrate-Resolved Spectral Fits to B 



Source 


Model Nh r 

(10 "cm ■') 


xS/dof 


Lx 
crgs/s 




Allowing PL norm and index to 


vary 




B (High) 
B (Med) 
B (Low) 


PHABS(POW) (5) l-5io2 

3 

i.U_Q 2 


0.84/36 


S.l+i'i X 10^4 
2.7+°-^ X 1034 

1.2iO-| X 1034 




Allowing PL norm and Nh to 


vary 




B (High) 
B (Med) 
B (Low) 


PHABS(POW) 5.0+f;-3 lA+'l^i 

O.M_Q 9 

22+_f 


0.83/37 


i X 1034 
2.4+0-5 ^ ^(,34 



Note. — Spectral fits to NGC 6652 B, split by countrate (see text). Errors are 
90% confidence for a single parameter. Luminosities are unabsorbed in erg s~^ for 
0.5-10 keV. 



